Estimation of qubit states in a factorizing basis 
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The optimal estimation of a quantum mechanical 2-state system (qubit) - with N identically 
prepared qubits available - is obtained by measuring all qubits simultaneously in an entangled basis. 
We report the experimental estimation of qubits using a succession of -/V measurements on individual 
qubits where the measurement basis is changed during the estimation procedure conditioned on the 
outcome of previous measurements (self- learning estimation). The performance of this adaptive 
algorithm is compared with other algorithms using measurements in a factorizing basis. 
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A question of fundamental and practical importance 
regarding the quantum mechanical description of the mi- 
croscopic world is: How can we obtain maximal infor- 
mation in order to characterize the state of a quantum 
system? Quantum states of various physical systems such 
as light fields, molecular wave packets, motional states of 
trapped ions and atomic beams have been determined 
experimentally with considerable precision 0. Acquir- 
ing complete knowledge about a quantum state would, 
of course, only be possible, if infinitely many copies of 
a quantum state were available and could be measured. 
More to the point, the initial question may be reformu- 
lated as the following task: Find a procedure consisting of 
a finite number of measurements yielding a state vector 
that best represents the (classical) knowledge possibly 
gained from any type of measurement of the quantum 
system under scrutiny. 

Determining an arbitrary state of a quantum mechan- 
ical two-state system (qubit) is of particular importance 
in the context of quantum information processing. In 
Rcf. H two identically prepared 2-state quantum sys- 
tems were considered with no nonlocal correlations and 
it was searched for the optimal measurement strategy 
to gain maximal information (difference of Shannon en- 
tropy) about this quantum state. It was strongly sug- 
gested that optimal information gain is achieved when a 
suitable measurement on both particles together is per- 
formed. Later it was proven that, indeed the optimal 
measurement for determining a quantum state - if two 
spin 1/2 particles are available - needs to be carried out 
on both particles together, i.e. the operator characteriz- 
ing the measurement does not factorize into components 
that act in the Hilbert spaces of individual particles only 
||. Moreover, an optimal estimate of the spin direction 
(the qubit state) of an ensemble of N identically prepared 
particles requires the application of such a nonfactoriz- 
able measurement operator. As a special case of optimal 
quantum state estimation of systems of arbitrary finite 
dimension the upper bound (N + 1)/(N + 2) for the mean 
fidelity of an estimate of N qubits was rederived in Ref. 

In particular, it was shown that finite positive op- 
erator valued measurements (POVMs) are sufficient for 
optimal state estimation. This result implied that an ex- 



perimental realization of such measurements is feasible, 
at least in principle. Subsequently, optimal POVMs were 
derived to determine the pure state of a qubit with the 
minimal number of projectors when up to N = 5 copies 
of the unknown state are available [^) . Still, the proposed 
optimal and minimal strategy requires the experimental 
implementation of rather intricate nonfactorizable oper- 
ators for a simultaneous measurement on all N qubits. 
First experimental steps towards entanglement-enhanced 
determination (N = 2) of a quantum state have been 
undertaken Q . Estimating a quantum state can also be 
viewed as the decoding procedure at the receiver end of 
a quantum channel necessary to recover quantum infor- 
mation (e.g. encoded as a unit vector) [0J|]. 

It was recently shown that quantum state estimation 
with fidelity close to the optimum is possible when a self- 
learning algorithm is used and measurements on N iden- 
tically prepared qubits are performed separately, even 
successively J9|. Here, we present, to our knowledge, the 
first experimental realization of a self-learning measure- 
ment on an individual quantum system in order to es- 
timate its state. The base of measurement is varied in 
real time during a sequence of N measurements condi- 
tioned on the results of previous measurements in this 
sequence. In addition, we compare the attainable ex- 
perimental fidelity of this adaptive strategy for quantum 
state estimation with strategies where the measurement 
base is either a predetermined one, or is randomly chosen 
during a sequence of N measurements. If a self-learning 
algorithm is employed to estimate a quantum state, then 
a suitable target function (here, the gain in the expected 
mean fidelity as described below) is maximized when pro- 
ceeding from measurement n — 1 to n. Under realistic 
experimental conditions possible errors have to be taken 
into account that may influence different measurement 
strategies differently. 

Here, the quantum mechanical two-state system un- 
der investigation is the S1/2 ground-state hyperfine dou- 
blet with total angular momentum F = 0, 1 of a single 
171 Yb + ion confined in a miniature Paul trap (diameter 
of 2 mm). The |0) = \F = 0) <-> \F = 1, m F = 0) = |1) 
transition with Bohr frequency ojq is driven by a quasires- 
onant microwave (mw) field with angular frequency near 
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lo = 2ir 12.6 GHz. The system is virtually free of deco- 
herence, i.e. transversal and longitudinal relaxation rates 
are negligible |l0|,[llj . Photon-counting resonance fluores- 
cence on the Sx/2(F=1) <-> P 1/2(^=0) transition driven 
by a frequency-doubled Ti:sapphire laser at 369 nm 
serves for state selective detection. Optical pumping into 
the \F = l,mj? = ±1) levels during a detection period is 
avoided when the E vector of the linearly polarized light 
subtends 45° with the direction of the applied dc mag- 
netic field. The light is detuned to the red side of the 
resonance line by some 20MHz in order to laser-cool the 
ion. Optical pumping the ion into the metastable 2 D 3 / 2 
level is prevented by illumination with light at 935 nm of 
a diode laser that retrieves the ion to the ground state via 
the \D 3/2 , F=l) |[l/2]i /2 , F = 0) excitation. Cool- 
ing is achieved by simultaneously irradiating the ion for 
100 ms with light from both laser sources and with mi- 
crowave radiation. This is done before each succession of 
measurements that consists of preparing and measuring 
a qubit state N times. 

In the reference frame rotating with lv, after apply- 
ing the rotating wave approximation, the time evolu- 
tion operator determining the evolution of the qubit ex- 
posed to linearly polarized mw radiation reads U(t) — 
exp [— |i (5a z + £la x )~\. The Rabi frequency is denoted 
by Q and a z . x represent the usual Pauli matrices. Any 
pure state can be represented by a unit vector in 3D 
configuration space (Bloch vector): |0, <f>) — cos||0) + 
sin |e l< ^|l) and is prepared by driving the qubit with mw 
pulses with appropriately chosen detuning S = luq — w, 
intensity, and duration, and by allowing for free preces- 
sion for a prescribed time. Rabi frequency (fl = 3.47 x 2n 
kHz) and detuning (i5 = 107 x 2ir Hz) of the mw radiation 
are determined by recording Rabi oscillations over 4-8 pe- 
riods and by performing a Ramsey-type experiment with 
mw pulses separated in time. A measurement in a given 
direction is performed in two steps: First, a suitable uni- 
tary transformation of the qubit is performed effecting 
a rotation of the desired measurement axis onto the z- 
axis. Second, the qubit is irradiated for 2 ms with laser 
light resonant with the S!/ 2 (F=1) <-* P1/2 transition and 
scattered photons are detected, if state |1) is occupied. 

A self-learning measurement of the prepared 
qubit state consists of N sequences each com- 
prising i) the preparation of \9 prep , pr ep), h) per- 
forming a projective measurement in the basis 
(1 0m, 0m) 11 ,\6 m = w- m , (f) m = TT + m ) n ) , and iii) us- 
ing the result of this (n— l) th measurement to determine 
the basis of the subsequent n th measurement that max- 
imizes the gain of the expected mean fidelity ||. This 
third step will be detailed in what follows. 

After n — 1 sequences the density operator repre- 
senting the state to be estimated is given by g n -i = 
d9sm9 J Q 27r d(j) w n -x{6, 4>) \9, 0)(0, <p\. The normalized 
probability density distribution w n -i(0, 0) is updated af- 
ter each measurement using Bayes rule Q, i.e. if in se- 
quence n the system is measured in direction (0 m ,0 m ) 



the distribution is modified by the probability for this 
outcome 

,n nn , \ W«-l(0,0) | (fl m , m |fl, 0) | 2 

W n [0, <p\9 m , m ) = J- — , (1) 

Pn\t>m, 0m) 

where the probability p n (0 m ,0 m ) = (9 m , <f> m \g n -i \&m, 0m) 
to find the system in direction (# m , (f> m ) in the n— th mea- 
surement ensures correct normalization. 

The best estimate of the pure qubit state, 
\0est, 0est)n-i is obtained by maximizing the fidelity 

^n-l(0>0) = (0,<P\Qn-l\&,(f>), i-e. ^ n -l(0 CS t,0ost) = 

F°^. 1 = maxF„_i(S, <j>). In order to find the optimal 
measurement direction for sequence n, the expected mean 
fidelity after measurement n is maximized as a function 
of the measurement direction. Suppose in the n-th mea- 
surement the qubit is found in direction (0 m ,0 m ). Then 

F„(0,0|0 m) m ) - 

Cde'shxe' f 7r d0' Wn (0',0'|0 m ,0 m )|(0,0|e',0')| 2 , (2) 

Jo Jo 

where the expected distribution w n (6', (f>'\0 m , m ) is ob- 
tained from Bayes rule (eq. |l|). The optimal fidelity 
F° pt (9 m , 4> m ) is obtained by maximizing this function 
with respect to (9, <f). Measurement n is performed along 
a specific axis and the qubit might as well be found in the 
direction (0 m) 0m) Therefore, the expected mean fidelity 
after the n— th measurement is given by the optimized 
fidelities for each of the two possible outcomes, weighted 
with the estimated probability for that outcome: 

F n (9m, m ) = Pn {9 m , m )F„° pt (0 m , m ) 

+ P n (0m,0m)F° pt (0 m ,0 m ) . (3) 

The optimal measurement direction (0^,0^*) max i- 
mizes this function. 

The direction of the first (n — 1) measurement is of 
course arbitrary, since no a priori information on the 
state (u>o(0, 0) = jj) is available. The expected mean fi- 
delity in this case is F\ = 2/3, independent of (0 m , m )i- 
After the first measurement the symmetry of the proba- 
bility distribution wi(0,(j>) is reduced to rotational sym- 
metry around the first measurement axis. 

The expected mean fidelity now depends only on the 
relative angle a between the second and the first mea- 
surement direction and we find F2 — (1/2 + cos(a/2 — 
7r/4)/vl8). Thus, the optimal second measurement with 
a = tt/2 yields F° pt = (1/2+ l/\/18). After the second 
measurement W%(9, 0) is still symmetric with respect to a 
plane spanned by the first two measurements directions. 
Again, the optimal measurement direction axis is orthog- 
onal to both previous directions and we obtain F^ p = 
(1/2 + 1/VT2). The optimal directions of subsequent 
measurements (n > 3) do depend on the outcome of pre- 
vious measurements. For an estimation procedure com- 
prised of N sequences, we have calculated numerically 2 N 
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possible successions of directions {(6 m , 4> m ) n } and pro- 
grammed the computer interface that controls the exper- 
imental parameters to choose the optimum measurement 
direction online during an estimation procedure. Fig. [l] 
illustrates a succession of measurements that yield an es- 
timate of the initial state |# pr ep, 4> P rc P ) — K/4, 7r/4) em- 
ploying the self-learning algorithm. The probability den- 
sity w n (8, 4>) is shown on the surface of the Bloch sphere 
and the n th and optimized (eq. ^) (n + l) th measurement 
directions are indicated. 




FIG. 1. Probability distribution w n (9,cf>) on the Bloch 
sphere updated by Bayes rule, in a typical realization of 12 
successive measurements. Darker areas correspond to higher 
probabilities (scaled individually for each Bloch sphere); con- 
tour lines for w n — 0, 0.1, ... are given. The n th and the 
optimized (n + l) th measurement directions are indicated by 
the open and the solid arrows respectively. The white circle 
shows the prepared state \0 prcp , 4> P rc P ) ~ [t/4, n/4). 

The discussion so far is based on the assumption, that 
measurements are performed with perfect efficiency. This 
is obviously not true in a real experiment. In this para- 
graph we will discuss the influence of experimental imper- 
fections on the quality of state estimation. Since the Rabi 
frequency and detuning 5 are determined precisely with 
an error below 1%, the deviation of the prepared state 
and of the measurement axis from their anticipated di- 
rections is small and the resulting systematic error in the 
fidelity is negligible compared to the statistical error. If 
there were no background signal during a detection pe- 
riod, the observation of m > scattered photons in a 
single measurement would reveal the ion to be in state 
|1) with probability 1 — pi(0) close to unity. (The prob- 
ability pi(m) to detect m photons follows a Poissonian 
distribution with mean value fh\ ps 5.) However, due to 
scattering off the ion trap electrodes and windows some 
photons will be detected even if the ion had been pre- 
pared in state |0) (also with a Poissonian distribution 
Po(m) with rhg ~ 0.2). In order to assign a given num- 
ber of photon counts in an individual measurement to 
the corresponding state of the ion, the threshold s is in- 
troduced: The probability rji to detect m > s photons, 
when photons are scattered off the ion (state |1)) is given 
by m = Em= s -Pi( m )- Analogously, n = Em~ioPo( TO ) 



for state |0). The functional relationship between ni and 
s is determined by the observed photon number distribu- 
tions pi(m). Since the detection efficiencies rji < 1, both 
a statistical and a systematic error are introduced into 
the measurements, as will be shown below. 

Using the average efficiency fj = (rjo + r]i)/2 and 
the efficiency difference An = (r/i — 7/o)/2, the prob- 
ability to find an "on" event (m > s) is given by 
P("on") = (2fj — l)Pi + (1 — fj) + An, and, analogously 
P("off") = {2fj-l)P + (l-fj)-Ar), where P = |(i|*)| 2 
and the |W) is the ion's state before irradiation with UV 
light. This effect of the measurement can be thought 
of as the distorting action of a quantum channel on 
the system's state followed by a perfect measurement: 
q — > (2fj — 1)q+ (1 — fj)I + An a z The channel acts as a 
depolarizing one characterized by the damping parameter 
1 — fj. The error introduced hereby is independent of the 
choice of the measurement basis and hence statistical. Ef- 
fectively the purity of the state (or equivalently the length 
of the Bloch vector \(a)\) decreases. The term in the fi- 
nal density matrix containing An systematically shifts 
the resulting state along the measurement direction. If 
an algorithm for state estimation is used that relies on 
measurements in fixed directions, for example in the x,- 
y- and z-direction, then the estimated state acquires a 
component parallel (or anti-parallel for A77 < 0) to the 
direction determined by the vector sum of the measure- 
ment directions. On the other hand, algorithms using 
measurement directions distributed over the whole Bloch 
sphere tend to cancel this error. This can be achieved 
with both the self-learning and the random algorithm. 
For experimental reasons we implemented only measure- 
ment directions on the upper hemisphere (i.e. 9 m < ir/2) 
and thus observe this systematic error for all algorithms 
if An 0. Choosing the threshold s op t such that An = 
eliminates this systematic, basis dependent error. When- 
ever an efficiency difference cannot be avoided, any algo- 
rithm can be made more robust against a systematic error 
in the state estimation by choosing measurement direc- 
tions such that their vector sum is close to zero. Note 
that An = (no bias in the estimation procedure) does 
not, in general, occur at that threshold that is required to 
make the most probable assignment to state |0) or |1) of 
a given number of detected photons. The threshold s' opt 
required for the latter would be at the intersection of the 
two distributions, yielding the maximum of the detection 
efficiency fj. However, the photon count distributions in 
our experiment yield s op t = s' opt = 2. 

We have studied the influence of the bias direction on 
the performance of all three algorithms. To this end A77 
was varied by changing the threshold s for the estima- 
tion of four different prepared states. Each state was es- 
timated several hundred times after 12 consecutive mea- 
surements for a given value of An. Fig. || shows that 
the dependence of the fidelity on An strongly varies for 
different states to be estimated. The curves in Fig. 2 in- 
tersect where the fidelity is independent of the prepared 
state. This intersection occurs at A77 = as is expected, 
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if the functional dependence of Arj on s is correct (de- 
termined independently using the experimental photon 
count distributions.) 



100 




FIG. 2. Mean fidelity attained with N = 12 successive 
measurements, optimized by the self-learning algorithm, as a 
function of the efficiency difference Ar/ for different prepared 
states (circles: j37r/4, 7r/4), diamonds: |37r/4, 37r/4), stars: 
Itt/4, 37r/4}, plus: |7r/4, 7r/4)). Each data point is averaged 
over 100-200 realizations. 



In addition to the fidelity optimizing adaptive algo- 
rithm, an "orthogonal" and a "random" one have been 
implemented for comparison. For the orthogonal algo- 
rithm equally many measurement are carried out in the 
x- y- and z- direction. The result of a succession of 4 mea- 
surements in each direction is evaluated by updating the 
probability distribution on the Bloch sphere using Bayes 
rule as is done for the adaptive algorithm. The random 
algorithm is realized by employing 2 N randomly gener- 
ated directions instead of the 2 N optimized directions as 
described above for the self- learning measurement. 

Table | shows the fidelities for these three algorithms 
together with the respective values expected from the- 
ory. The experimental fidelities for each algorithm are 
obtained at the intersection of four curves at = 
corresponding to the estimation of four different initial 
states as described above. The attainable fidelity is lim- 
ited by experimental imperfections, i.e. by the finite de- 
tection efficiency fj = 97%, and most notably by the im- 
pure preparation of state |0) at the beginning of each 
sequence of measurements (%> r0 p = 89%). The orthogo- 
nal algorithm in addition suffered from light-induced de- 
coherence PJ. These imperfections reduce the purity of 
the state, i. e. the length of the Bloch vector | (a) | . This is 
accounted for in the theoretical values (table , that are 
average values obtained from numerically simulating each 
algorithm 10.000 times. It should be emphasized that the 
fidelities given are average values valid for measurement 
sequences with N = 12. If instead, the information of all 
sequences (typically N x 100 = 1200) is used for state 
estimation, then the experimental fidelity is better than 
99%. 



Algorithm 


( J F>Ex P . 


K*>l 


(JWG <<?)|) 


self-learning 


85.0 ± 0.6 


74.8 ± 2.1 


85.4 ±0.7 


random 


81.9 ±0.6 


73.4 ±2.1 


81.9 ±0.7 


orthogonal 


67.8 ±1.1 


40.4 ±4.6 


70.2 ±2.1 


orthogonal 




74.1 


83.6 



TABLE I. Experimental mean fidelities (taken from the 
data plotted in Fig. |i| at: At; = 0), overall length of Bloch 
vector | (a) | and mean fidelities expected from theory, when 
the length of the Bloch vector is taken into account. For 
comparison the expected value for the orthogonal algorithm 
is also given for a |(<r)| as realized with the other algorithms. 



The method and results presented are not restricted 
to a particular realization of qubits. The identification 
of imperfections in our experiment show that the fideli- 
ties obtained are currently limited mainly by the impure 
preparation. Here, a significant improvement seems fea- 



sible: 



?7pre 



> 



99% would lead to mean fidelities (with 



N = 12 qubits) better than 90%, close to the upper 
bound of 93%, attainable with an entangled measure- 
ment. 
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